TUM/T39-97-6 
UFTP 439/1996 

Perturbative Part of the Non-Singlet Structure Function F 2 in the Large-iVV 

LimitQ 



L. Mankiewic^ 
Institute for Theoretical Physics 
TU-Miinchen, D-85747 Garching, Germany 

and 

M. Maul, E. Stem 
Institute for Theoretical Physics, Johann Wolfgang Goethe- Universitat, 
D- 6 '90 '54 Frankfurt am Main, Germany 



We have calculated MS Wilson coefficients and anomalous dimensions for the 
non-singlet part of the structure function F 2 in the large- Np limit. Our result 
agrees with exact two and three loop calculations and gives the leading Np 
dependence of the perturbative non-singlet Wilson coefficients to all orders 
in as- 



PACS: 12.38. Cy, 13.60.Hb, 12.38.-t 

Keywords: deep inelastic scattering, structure function F 2 perturbative coefficients, 
anomalous dimensions 



1 Work supported in part by BMBF 

2 On leave of absence from N. Copernicus Astronomical Center, Polish Academy of Science, ul. 
Bartycka 18, PL-00-716 Warsaw (Poland) 



1 



High precision of current experimental data on deep-inelastic lepton-nucleon scatter- 
ing makes it necessary to extend perturbative calculations to higher and higher orders. 
Although in principle QCD perturbation theory is well understood, in practice its tech- 
nical complexity has limited the existing state-of-the-art next to next to leading order 
(NNL) calculations only to low Mellin moments M^, N = 2,4,6,8 of the structure 
functions F 2 and F L 0, g|. 

Recently, it has been realized that in the large- Np limit of QCD corresponding cal- 
culations can be performed exactly to all orders || ||. Although it is perhaps not easy 
to argue about the importance of results obtained in this limit for the real world, they 
are certainly useful either as an independent check of complicated computer algebra 
algorithms used in the exact calculations, or as a starting point for further approxi- 
mation procedures, such as 'Naive Nonabelianization' (NNA) j| or 'Asymptotic Pade 
Approximation' |J. 

The goal of the present paper is to analyze moments M-2,n{Q 2 ) of the flavor non- 
singlet (NS) component of twist- 2 part of structure function F 2 (x,Q 2 ) 

M 2 . N (Q 2 )= f dxx N ~ 2 F 2 NS (x,Q 2 ) iV = 2,4,..., (1) 
Jo 

defined through the well known decomposition |7[| of the hadronic scattering tensor of 
unpolarized deep inelastic lepton nucleon scattering in terms of the structure functions 
F 2 (x,Q 2 ) and F L (x,Q 2 ). 

W^p,q) = ±- Jd 4 ze^(p\.Uz),UO)\p)= U v - ^ ±-F L (x,Q 2 ) 



V (p ■ q) v ■ q j 2x 



(2) 



According to the Operator Product Expansion || one can arrange A4 2} n(Q 2 ) as a prod- 
uct 

M 2 , N (Q 2 ) = C 2 , N {a s {Q 2 ))A N {Q 2 ) (3) 

of the Wilson coefficient C 2j n(c^s(Q 2 )) an d the spin-averaged matrix element A N {Q 2 ) of 
a spin-N, twist- 2 operator 

{pl^^HD^ . . . iD^ N} i/j\p) p ~ n = p {w-pwv> A N (/j, 2 ) . (4) 

In the above equation {/i . . . v} indicates twist-2 projection i.e., symmetric and traceless 
combination. Validity of equation (§) requires that both renormalization and factoriza- 
tion scales have been set equal to the virtuality of the external photon Q 2 . In addition, 
it has been implicitly understood that squares of quark charges are properly included on 
the right-hand side of Eq.(^j), and that flavors of quark-operators if) have been combined 
to yield the difference between u- and d-quark matrix elements. 
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The Q 2 -dependence of matrix elements An(Q 2 ), written as a solution of the renor- 
malisation group equation, 

A N {Q 2 ) = A N {f) • exp [ r m das' ^^'n . (5) 

is determined by a set of anomalous dimensions 7at(o:s'(<3 2 )). The perturbative expansion 
of C2,n and 7^ results in 

00 

c 2 , N (as(Q 2 )) = i + E^V, +1 , 

3=0 

CO 

lN(a s (Q 2 )) = Y.1n< +1 , (6) 

3=0 

where we have introduced a shorthand notation a s = as< p \ In the following we present 
results for the Wilson coefficients C2,n(®s(Q 2 )) an d anomalous dimensions 7jv( a s(Q 2 )) 
obtained in the large- Np limit, Np — > 00, keeping a s ■ Np = const. 

Technically, we have found it advantageous to follow the approach developed in J|, ||] 
and calculate the generating function G(u; N) 

2 », V d j 
1 N *) ** 



Giu-^ + OiNiT 1 ), (7) 

for the Wilson coefficients C^jv i n the large- Np limit. As explained in details e.g. in |||, f|], 
the dominant contribution in this limit arises from an arbitrary number of fermionic loop 
insertions into the gluon propagator. As a consequence, the problem can be reduced to 
calculation of the first-order radiative corrections using modified Landau-gauge gluon 
propagator 

D ^ k > u ) = ( _ fc 2)i+u - ipr) • ( 8 ) 

In calculations of gauge invariant quantities, such as in the present case, the longitudinal 
part does not contribute and can be dropped. The only complication arises due to 
renormalization - the 'bare' function Gb{u; N) acquires a singularity at u = and has 
to be renormalized by an appropriate counterterm. As discussed in in minimal 

subtraction schemes such as MS the counterterm can be obtained in a compact form. 
As a result the renormalized generating function Gr(u; N) can be written as 

G R (u;N) = G B (u;N) + G (u;N)/u , (9) 

where Gq(u; N) is determined in terms of expansion coefficients gk{N) of another function 
Go(u; N) in the variable u around the origin 

00 

G (u;N) = Y,9k(N)u k , 

fc=0 

00 1 

G (u;N) = Y.^9k{N)u k . (10) 

fc=0 
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The merit of this complicated and apparently indirect construction lies in the fact that 
the function Gq(u; N) can be found directly from the calculation of one-loop radia- 
tive corrections computed in d-dimensions using the modified propagator (|8]), while the 
anomalous dimension follows from the corresponding counterterm in a standard way. 
An interested reader should consult Refs. [|], |9|, in particular Appendix A of Ref. j|, 
for details. To compute G R (u; N) and 77V we have calculated the Feynman diagrams 
depicted on Figure [I], arriving at 

G R (u;N) = G^u-N)+G B R {u-N) + Gl(u-N) + G D R (u-N) 



G%(u]N) = G B R (u-N) = 2C F [exp(-C)] u - 1(1 



G c R {u-N) = 12C F [exp(-C)] u 



G c Q {u-N) 



uT(3-u)(l + u) 
-2C F [eM-C)} T{3 _ uni + u)mKu + N) 

f 2 r rexnf nv T{1 ~ u) T{u + k) 

-^2C F (eM-C}) T{3 _ uni + u)m{u + k) 

+ g ^(ex P [-cir tir(3 _ u)r(1+u)r(jfc)(M+jfc+1) 

, G$(u; N) 
u 

N 



G^N) = r < 4 + 2 "» 



3 r(i- M )r 2 (2 + M )r(i + u)(u + fc) 

C F r(4 + 2w) 

Tr(i -w)r 2 (i + M )r(3 + M) 

r(l - u)T(N + u) 



-2C F [exp(-C)} u 



r(3 -u)(l + u + N)T(N)T(l + u) 

r(i - m)(m 2 - 3m + 2)r 2 (iv + u ) 



«r(3 - u)r(l + u)T(N + u + 2)T(N) 

GC(u;N) 



u 

C F r(4 + 2n)r(iV + n) 

~3~r(l - u)T 3 (l + u)T(N + u + 2) 



G R (u; N) = -2C F [eM-C)T ^ " ^ + " 



«r(3 - u)r(i + u)r(iv) 



3 r(i — u)T z {i + -u)r(3 + -u) 
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for the generating function Gn(u;n). In the above formula Cp = 4/3, and C denotes 
the finite part of the quark loop insertion into the gluon propagator, Cj^ = —5/3. The 
relation between G (u; N) and G Q (u; N) is given by equation (|T0l). For the case of the 
anomalous dimension 7Ar(as) our result reads 



7nM 



7nM 



b i x r 1 r(4 + 2s) 

lN (a s ) - -^- r(1 _ s)r2(1 + s)r(3 + s) 

fi r( 4 + 2g ) 

S F ^ 2 3T(l-s)T 2 (2 + s)T(l + s)(s + k) 
1 T(4 + 2s)T(N + s) 

as F 3r(i-s)r 3 (i + s )r(iv + s + 2) 

1 T(A + 2s) 

5 F 3r(i- s )r 2 (i + s )r(3 + s) ' 



(12) 



where s = ^Np^j a s . Terms labeled as A,B,C and D in Eqs. (|TT|) and fll2"D corre- 
spond to contributions of graphs depicted in Figure |l|. We quote results obtained in 
the Feynman gauge i.e. contributions from the longitudinal part of the propagator (|8|), 
which cancel in the sum of all graphs, have been neglected. In the case of the anomalous 
dimension 7Ar(a s ) one can put together various terms in ( |T2| ) and obtain a more compact 
expression 



lN{a s ) 



dsCw— 



1 T(4 + 2s) 



3r(i-s)r(i + s ) 3 

1 1 2 » I 

{s + 2)(s + 1) (s + N + l)(s + N) + (1 + s) 2 ^ 2 fc + s 



(13) 



(— a s , which agrees with result obtained earlier by Gracey (T^ . 



Now, expanding the right-hand side of (|T3|) in 

7Ar(a s ) = feiCFCts + C F 



(14) 



one can find explicit results for the anomalous dimensions in any fixed order of the per- 
turbative expansion in the large- iVV-limit, see Ref. [[[(J. Similarly, using (]?]) one can 
derive from Eq. ([TI]) fixed order predictions for the Wilson coefficients for the structure 
function F 2 in this limit. Note that large- Np generating function for the Wilson coeffi- 
cients for the non-singlet part of the structure function Fl(x,Q 2 ) is known from Refs. 

hi m. 
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In Ref. Wilson coefficients for Fl and F2 have been calculated for the even 
moments N = 2 . . .8 exactly up to the third order in a s . The leading Np terms at the 
third and the fourth orders for an arbitrary iV-th moment read 



^(3) 
°2,iV 



r»(4) 

°2,JV 



c F --m 



3 

10 1 

T (1 + iv) 4 



9517 _ J_ 5659 + 108 C(3) 

216 ^ U 81 1 



1 



1 

+ 8l 



-2041 + 108 C(3) 



+ (gC(3) + 

392 1 



4357 



9 1 + N 
, 19 2 
+ ("7T + 



162 A^ 3 
)Si(iV-l) 



+ 



N 
32 1 



30 

AT2 



N 
134 



9 



164 1 



9 



iV 
2 



+ 



368 

IT (1 + N) 2 
10 1 . 4955 1 
r iV 
56 1 
"3 



1 



146 

IT (1 + 7V)3 



iV 3 3 iV 4 + 



+ iV I \ 



+( 



N 2 
1 2 1 

3 ~ 3 iV 
535 2 



16 1 

T N 
2 1 



+ 



28 



1 



166 



31 + F 
16 1 

N 2+ liN + (1 + N) 2 
4 1 4 1 



(1 + N) 2 
\S 1 (N-1) 3 

2 28 
3 



3 1 + N 



(1 + N) 3 3 (1 + N) 2 
Si(JV - 1) S 2 (iV - 1) 

)Sx(iV-l) 2 



1 + N 



20 



v 9 

N-l 

+E 

k=l 

N-l 

+E 

fe=i 

JV-l 

+E 



, S 3 (iV- 1) S 4 (JV- 1) 

3 N 3 1+JV' 3V J 3 K J 

2 (-146 fc 4 - 139 k 3 - 8 k 2 + 3 k - 18) Si(fc) N 

9 (A; + 1) 3 P , 

2 (18 A; 4 + 36 k 3 + 18 fc 2 ) Si(Jfe) S 2 (Jfe)\ 



V 9 

2 (16/c 2 + 7/c 



(k + l) 3 k 2 
3)S^k) 2 ^4^) 



Mfc + 1) 2 
^2 (16 /fc 2 + 7 A; - 3) S 2 (fc) 



k=l 



--Np 
3 



k (k + l) 2 

386 1 



^E 

fc=i 

AT-l 

+ E 

k=l 

325 



3 jfe + 1 

8 S 3 (fc) 
3 k + 1 

1 



+ 0(N F ) 



1_27C(3) 
27 
1 



3 (1 + AT) 3 
2068 1 5713 
27 



AT2 

452200 + 5940 C(3) 



5V-1 

6 (1+iv) 4 1 v£ 

27C(3) 4955 1 

7V2 ' 



-6 



S 4 (fc). 



(1 + N) 2 
2430 C (4) 



k + rj 
83 265 
N 3 



jCO) 



1 

iV 4 



34883 1 



27 iV 2 N 3 6 N 4 162 
1 540 C(3) + 154360 + 2430 C(4) 



1620 

23 1 23 

+ y + y 

2077 1 

iV 



1 + Af 
+ 



1620 



+ 



27 
19 
Y 



(1 + AO 5 
12 

+ 



N 

A/ s 10 A/ x 25279 12 32 164 
-3 C(4) + - C(3) + + — - — + 



+ ~C(4) 



56 



3 
392 



324 ' N 4 
5731 



N 3 



1 

N 2 



(1 + NY 



6 16 

N 2 + N + (1 + N) 2 



(1 + iV) 3 
6 28 



+ 



1 



1 + N 



3 (1 + A0 2 ' 27 1 + N 
Si(JV-l)S 2 (N-l) 



)S X (JV-1) 



6 



+ (-i-^ + ra' s ^- 1 ' 4+ (-I-^ + ra) s ^- 1 ' 2 

203 6 16 82 1 6 28 196 1 2 

+( ~T2" ~ iv 3 + n 2 ~ T iv + (i + iv) 3 ~~ (i + at) 2 + T" TTiv^ ( ~ ^ 

,19 2 16 1 2 28 1 3 

+( ""6~ + W 2 ~ T n - (TTivF + T T+iv } l( " } 

+(^ + | - Si(JV - l) 2 S 2 (N - 1) + (-2 - 1 + ^) Sx(JV - 1) S 3 (iV - 1) 

, 10649 6 16 82 1 6 28 

+(— T^T " 2 C 3 + ^- ^ + ^T77- 



108 A'" 3 N 2 3 N (1 + iV) 3 (1 + iV) 2 

)S 2 (iV-l) + (— + — ---- - — + - — -) S 3 (iV - 1) 



3 1 + JV' v 3 N 2 3 N (1 + iV) 2 3 1 + iV 

^ (_\_ (146 fc 4 + 139 A: 3 + 8 k 2 - 3 fc + 18) S 2 (fc) \ ^ / S 2 (fc) 2 
+ k V 3 fc 2 (fc + l) 3 ) + k\ + 1 

+ ^ /l (146 A; 4 + 139 A; 3 + 8 A; 2 - 3 fc + 18) Si(fc) 2 \ 



fc=i 



A; 2 (fc + 1) ; 



^ / 1 (324 - 3128 fc 6 - 4929 fc 5 - 2427 fc 4 - 257 fc 3 - 333 fc 2 + 270 fc) Si (fc) ' 

^ /2 (16fc 2 + 7fc-3)S!(fc) 3 \ ^ Sx(fc) 4 _ 281971 
+ \3 fc(fc + l) 2 J + fr 1 A;+l 1728 

+ (-^ + ^- TT¥ )S 4 (iV-l) + 23S 5 (iV-l) 

^ /_ 1 (-216fc 6 - 648fc 5 - 648fc 4 -216fc 3 )Si(fc)S 3 (fc)^ 
k=i 

N-l 

+£ 

fc=i 

N-l 

+ £ 

fc=i 

JV-l 

+£ 

k=i 



27 (fc + i) 4 fc 3 ; 

'4 (16fc 2 + 7fc-3) S 3 (fc) \ 
3 fc(fc + l) 2 J 

' 1 (18fc 4 + 36fc 3 + 18fc 2 )Si(fc) 2 S 2 (fc) \ 
~3 (fc + l) 3 fc 2 J 

1 (864 fc 6 + 2106 fc 5 + 1458 fc 4 + 54 fc 3 - 162 fc 2 ) Si(fc) S 2 (fc) ' 
~27 (fc + l) 4 fc 3 



+ 0{N 2 F ) , 

(15) 



where S a (N) = I]^L 1 (l/fc a ). The corresponding results for the longitudinal structure 
function F L read 



L,N = \-o ly P 



1 + N 



18 3(1 + N) (1 + iV) 2 V 3 ( 1 + N ), 



+Sf(N - 1) - S 2 (N - 1) 



+ 0(7V F ) , 



7 



c 



(4) 
L,N 




+4Sf(N - 1; 

+8S 3 (N-1) 



(l + N) (1 + N)\ 
12 



38 



'1 




S\{N 



+ 12Sx(N 



C 



(5) 
L,N 



+ 



+ 



+ 



--N F 
3 



/19820 
v 27 

152 

812 

IT 

3 
812 

IT 

/304 

VlT 



c 



+ 0{Nl) , 
69766 



l + N 



81 



1624 



+ 



3(1 + N) 
18 )5 2 (iV 



19820 
27(1 + AT) 
304 

- 96- 



+ 



1624 



304 



3(1 + iV) 2 
1 



;i + ny 



+ 96 



(l + NY 



[l + N) 

152 
l + N 
16 



1 + iV) 2 (1 + A^) 3 
- 1) + 32S 3 (N - 1 



+ 



48 



(1+iV) 5 



245 2 (iV- i; 



Si(N- 
Sl(N 



+ ( ] f- Y f^)s!(N-l)+ASt(N-l) 



l + N 
152 

(TTiv) 

32 
1 + iV 



+ 



48 



1 + AT) 2 
5 3 (A^ - 1) - 245 4 (iV - 1) 



S 2 (N - 1) + 12S 2 2 (N - 1) 



+ 0(iV|) 



(16) 



respectively. We have checked that the first, second (not quoted here) and the third 
order coefficients agree with the leading- Np terms extracted from Ref. 0. Results for 
higher order terms are too long to be written down explicitly. Instead, we have prepared 
a Maple program for the numerical evaluation of Wilson coefficients in the large-TV^ 
limit, available on request, which allows to compute them using results of [10, 11] and 
of the present paper. 

Finally, in Tables [I] and |2] we have compared the NNA approximants for the coefficient 
functions C2,jv(a s ) with the exact results obtained in Ref. M, and evaluate the NNA 
prediction for the 0(af) terms. It is seen that although typically the NNA procedure 
predicts correctly the magnitude of the perturbative coefficients, with increasing preci- 
sion as A^ becomes larger, the numerical accuracy is not optimal. We have also checked 
that a similar procedure applied to the anomalous dimension 7Ar(a s ) gives much worse 
results which is probably connected to the fact that in the latter case the perturbative 
expansion is not dominated by renormalons. 
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Nf ~- 


= 3 


Exact results P 


NNA approximants 


N = 


2 


1.69377 a 2 + 1.42209 


71.99999 a 2 s + 1099.02 a 3 


+ 26193.11796 a* 


N = 


4 


91.3797 + 1675.76 


229.3368 a 2 + 4256.46 a 3 s 


+ 103655.6199 a* 


N = 


6 


218.3596 a 2 + 5004.63 


378.1762 a;; + 7775.28 


+ 200938.7405 a 4 s 


N = 


8 


357.0330 a 2 + 9357.69 


511.9851 a 2 + 11283.4 


+ 306602.1048 a 4 s 



Table 1: Comparison of the NNA approximants to the exact results of the coefficient function 
C2,n(cl s ) obtained in @ up to order 0{a 3 ) for N F = 3. The last column contains also the 
NNA prediction for the 0(a 4 s ) terms. The 0(a s ) corrections agree of course exactly. 



Nf = 


= 4 


Exact results 


NNA approximants 


N = 


2 


-3.63957 a 2 - 169.747. a 3 s 


66.66666 a 2 + 942.230 a 3 s 


+ 20792.94153 a\ 


N = 


= 4 


74.3918 a 2 + 901.570 a 3 


212.3489 a 2 s + 3649.22 a 3 


+ 82285.17302 a A s 


N = 


6 


190.3465 a 2 s + 3454.66 a 3 s 


350.1631 a 2 s + 6666.05 a 3 


+ 159511.6507 a 4 s 


N = 


8 


319.1081 a 2 s + 6973.59 a 3 s 


474.0603 a 2 s + 9673.72 a 3 s 


+ 243390.6360 a 4 s 



Table 2: Comparison of the NNA approximants to the exact results of the coefficient func 
tion C 2) Ar(a s ) obtained in up to order 0(a 3 ) for N F = 4. The last column contains also 
the prediction for the 0(a 4 s ) terms. The 0(a s ) corrections agree of course exactly. 
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— + -o ♦ oo ♦ ooo ... 

Figure 1: Graphs A, B, C and D which contribute to the calculation of the perturbative 
part of F 2 (x,Q 2 ) in the large- N F limit. 
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